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We present here analytical solutions of General Relativity that describe evolving wormholes 
with a non-constant redshift function. We show that the matter that threads these wormholes is 
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'. I. INTRODUCTION 

o . 

In the last two decades, there has been a revival on the study of classical wormhole solutions either in General 
Relativity (GR) or in alternative theories of gravitation. A cornerstone in this revival has been the paper of Morris 
and Thorne |l[ in which, they derived the properties that a space-time must have in order to hold up such a geometry. 
Let us recall that the most salient feature of these space-times is that an embedding of one of their space-like sections 
in euclidean space displays two asymptotically flat regions joined by a throat. There are several reasons that support 
£N| , the interest in these solutions. One of them is the possibility of constructing time machines |2|-|[ which violate 
Hawking's chronology protection conjecture || . Another reason is related to the nature of the matter that generates 
this solution: it must be "exotic", i.e. its energy density takes negative values in some reference systems. This 
entails violations of the weak energy condition (WEC) which, although not based in any strong physical evidence 
||, is deep enough in the body of relativists' belief so as to abandon it without serious reasons. These violations 
are in fact encoded in the evolution of the expansion scalar, governed by Raychaudhuri's equation [Q. As WEC is 
^jH, violated for any static spherically symmetric wormhole and has recently been shown to be the case also for any static 
wormhole in GR ||, the attemps to get around WEC violations have led two different branches of research. Firstly, 
there have been an increasing number of works in non-standard gravity theories. Brans-Dicke gravitysupports static 
wormholes both in vacuum (|] and with matter content that do not violate the WEC by itself [Jl(|. Also, there 
. -h [ are analysis in several others alternative gravitational theories, such as R + R 2 |pd| , Einstein-Gauss-Bonnet Jl2| and 
Einstein-Cartan models |fl3|| . In all cases, the violation of WEC is a necessary condition for static wormhole to exist, 
although these changes in the gravitational action allows one to have normal matter while relegating the exoticity 
to non-standard fields. This lead to explore the issue of WEC violation in non-static situations and see which is the 
case for dynamic, time-dependent wormholes. The first of this kind of analysis, due to Roman fbl|j , was devoted to 
the case of inflating Lorentzian wormholes. These are wormholes of the Morris-Thorne type embedded in an inflating 
background, with all non-temporal components of the metric tensor multiplied by a factor of the form e 2x ', being \ 
related to the cosmological constant. The aim of Roman was to show that a microscopical wormhole could be enlarged 
by inflationary processes. However, this construction also violate the WEC. The other works, due to Kar and Kar and 
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Sahdev 15 16|], pointed at testing whether, within classical GR, a class of non-static, non- WEC- violating wormholes 
could exists. They used a metric with a conformal time-dependent factor, whose spacelike sections are R x S 2 with 
a wormhole metric, like the ones analyzed in [Q. In this work, we shall extend this last analysis by studying more 
general cases of evolving wormholes, given by a generic Morris-Thorne metric, in the presence of matter described by 
a given stress-energy tensor. To proceed further we shall make two ansatze. Firstly, we give a specification of the 
redshift function $ of Morris-Thorne consistent, for any hypersurface of constant t, with asymptotic flatness and a 
metric free of horizons, as was the one made for the static case in pJV- Then we shall show how to choose a particular 
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functional form for the trace of the stress-energy tensor matter in order to get different solutions for the conformal 
factor. At this point, we shall derive an analytical solution which we use later to analyze the WEC violation scenario 
and some human tranversability criteria. 



II. GEOMETRICAL AND TOPOLOGICAL FEATURES 



We adopt the following diagonal metric for the space-time: 



ds 2 



n( t ) {-e^dt 2 + e 2A ^dr 2 + r 2 dfi 2 } 



(1) 



with O(t) a conformal factor, finite and positive defined throughout the domain of t, and <iil|, the S 2 line element 
as usual. When Sl(t) is constant, this general metric stands for the ones firstly analyzed in the work of Morris and 
Thorne JT|. When e 2 * = 1, it reduces to the particular case studied by Kar p5|j . In the spirit of |jl7), we make the 
Ansatz $ = —a/r, where a is a positive constant to be determined. This choice guarantees that the redshift function 
is finite everywhere, and consequently there is no event horizon. For the stress-energy tensor we take that of an 
imperfect fluid [] f|) , 



Tuv = {p + P) t^£ v + P 9fiu + Ap 



2 \9iw 



(2) 



where p stands for the energy density, p for the isotropic fluid pressure, and q for the energy flux. The anisotropy 
pressure Ap is the difference between the local radial (p r ) and lateral (p±) stresses. With the foregoing assumptions, 
the nonvanishing Einstein's equations are 
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In the expressions given above, dashes denote derivatives with respect to r while dots, derivatives with respect to t. 
It seems convenient to assume that the trace of the stress energy tensor will be a separated function of t and r, 
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After separating variables, equation (fa) can be rewritten as, 



n 3 in 

'fi + 2 I fi 



- T = C 



(8) 



1 In what follows, greek indices run from to 3, parenthesis denotes symmetrization, £ M = \f—g~oo6° represents time-like unit 
vectors and f M = ■JgnS* space-like unit vectors in the radial direction respectively. 
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where C is an arbitrary constant which for simplicity, we shall take equal to 0. Performing the change of variables 
A = lnz, equation (|^) transforms into a Bernoulli equation, 
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The additional change y = z~ 2 = exp [— 2A], allows for the general solution to be written as 

_ 2A _ <j> 2 + 9</> 3 + 32</> 5 + ( 57/2 + e 2 ^ /C )0 4 
V ~ e ~ (20 + 1)5 ■ l 11 ) 

with <p = r/a, and K, an integration constant to be determined. Note that, far from the wormhole throat, 

lim e- 2A = 1, (12) 

r — >oo 

or equivalently, 

dz 

lim — = 0, (13) 

i — >oo dr 

where z(r) is the embedding function p]]. Moreover, for any hypersurface of constant t the mouths of the wormhole 
need to connect two asymptotically flat spacetimes, thus the geometry at the wormhole's throat is severely constrained. 
Indeed, the definition of the throat (minimum wormhole radius), entails for a vertical slope of the embedding surface, 

lim ^= lim ±y/e 2A - 1 = oo; (14) 
th th 

besides, the solution must satisfy the flaring out condition Q, which stated mathematically reads, 

A'e- 2A 

- (i_ e -2A )2 >0 - ( 15 ) 
Equations ( p^ ) and (|l^) will be satisfied if and only if 

lim er 2A = 0+, (16) 

th 

thus, in order to fix the constant /C, we we must select a value for the dimensionless radius of the throat ((f> th > 0) such 
that Eq. ([l6]) be satisfied. Nevertheless, the absolute size of the throat also depends on a. As an example let us impose 
4> th = 2 which, using (p"T|), sets K. — — 97.25 e _1 and r th —2 a. The aforementioned properties of A, together with the 
definition of (f> and f2, entail that the metric tensor describes, for il(t) = constant, two asymptotically flat spacetimes 
joined by a throat. Thus, equation (|ll]) represents an analytical solution for an evolving wormhole geometry, which, 
due to the choice we made for the trace, is independent of fl(t). To obtain the complete behavior of the metric, we 
have to define T(t) in order to solve for f2(i) in equation (||). We take instead Eq. (|§|) as the definition of T(t) in 
what follows. 



III. AVOIDING WEC VIOLATIONS 



Having at our disposal an explicit wormhole solution, given in the form of y (</>), we can now turn to the question of 
whether it entails WEC violations. Of course, we have to note that, in the general case, it will ultimately depend on 
the specific choice of £l(t), but we use here the fact that, due to the form in which the solution for y((f>) was computed, 
it is independent of the explicit form of the conformal factor. Thus, we are able to study different functional forms of 
ri(t) for the same y(<j>). The WEC asserts that for any time-like vector v^ 1 : 



T^v v v» > 0. 



(17) 



The physical significance of this condition is that the observed energy density, as measured by any time- like observer, 
has to be positive. The connection of WEC violations with the existence of wormholes geometries has been widely 
studied before |,||Jl|. In the case of evolving wormholes, Kar Jl5| has demonstrated that there exist wormhole 
solutions which do not violate the WEC for a specific choice of the redshift function. His analysis is indeed general, in 
the sense that even without having worked out any explicit wormhole solution, he could conclude that there exist some 
Lorentzian evolving wormhole geometries with the required matter not violating the WEC. We must recall however 
that Kar dealt with the particular case <& = 0. We pursue here to refine his analysis, by means of the solution for the 
metric component gn given in the previous section. With the stress-energy tensor of equation (^), WEC entails 

P>0 p+p r + 2q>0 p + p±>0. (18) 

Starting from the Einstein equations, the three previous inequalities can be written as follows: 

A(t)+9i(r)>0 V(t,r) (19) 

f2(t)+g 2 (r) + 2q(r,t)>0 V(t,r) (20) 

fa(t)+gs(r) >0 V(t,r) (21) 

where the functions fi and gi are as in Table |[ It is a trivial exercise to write all radial functions in terms of the 
variable by means of the definition <j> — r/a. After this is done, it is seen that these functions display an explicit 
dependence on the a parameter and on </> th , through K-. The radius of the throat will be fixed when a and <p th are 
simultaneously known. These radial functions are plotted in Figs. 1,2 and 3 for different values of a and for </> th = 2. 
Note that, at or near the throat, some of them are negative. This just reflects the fact that no static wormholes 
can be built, within GR, without using exotic matter. Far enough from the throat, all functions tend to zero, being 
strong the dependence on a of the shape of the curves. Also from the figures, we may note that if we demand that 
the time-dependent functions f\ , fi and be greater than the modulus of the minimum of each of the corresponding 
radial functions, all WEC inequalities will be fulfilled. In the case of the second inequality, the term coming from the 
energy flux must be taken into account to determine the non-exotic region. In Table ||, we provide several examples of 
possible solutions for f2(i). The procedure to obtain these particular functions as solutions for the metric is, as stated, 
to define a suitable functional form for the temporal part of the trace of the stress-energy tensor for matter, i.e. T(t). 
This values are also given together with the values of the temporal functions fi for each choice. This entries would 
allow comparison with WEC inequalities in order to search for constraints on the values of the constants involved 
in each choice of to be fulfilled in order to obtain non-exotic matter at the throat of the evolving wormhole. 

Finally, the domain of t for which fl is finite and without zeros is also given in Table |l[ The comments made by 
Kar pl| when he considered different Q(t) functions in his scheme are totally applicable here. For instance, we also 
have exponentially expanding or contracting phases -for the first two solutions- and a wormhole -like universe (with 
a bang and a crunch) for f2(t) = sin 2 (cut). In our scheme, the particular choice of Sl(t) will determine the energy flux, 
through the field equation (^). The sign of the energy flux (which depends on Cl) will decide in turn if the wormhole 
is attractive or repulsive, in the sense explained in O], For large enough and positive q, the second WEC inequality 
will be fulfilled easily. Morover, this would imply a quick process of expansion which in turn favors a possible trip by 
diminishing the tidal forces, as we shall show in the next section. 

IV. TRAVERS ABILITY CRITERIA 

We shall deal now with the possibility that an explorer might enter into the kind of evolving wormhole studied 
here, pass through the tunnel, and exit into an external space again without being crushed during the journey. We 
are not going to discuss here the stability of the wormhole, being this beyond the scope of this paper. Instead, we 
content ourselves with an analysis of the tidal gravitational forces that an infalling radial observer must bear during 
the trip. The extension to a non-radial motion follows from the recipe given in |l9) (Chap. 13). We recall the desired 
traversability criteria concerning the tidal forces that an observer should feel during the journey: they must not 
exceed the ones due to Earth gravity [ pp[ . From now on, the algebra will be simplified by switching to an orthonormal 
reference frame in which g^t, = rj^p. In terms of this basis, the proper reference frame is given by, 

e 6' = 7 e t - =F 7 e? ej, = =p 7 e? + 7 f3 e~ t e 2 , = e § e§, = e; (22) 
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where 7 and /3 are the usual parameters of the Lorentz transformation formula. Thus, stated mathematically the 
traversability criteria read, 
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(24) 



Note that, contrary to the static case, as the geometry is changing with time, the tidal forces that the hypothetical 
observer feels also depend on time. The first of these inequalities could be interpreted as a constraint on the gradient 
of the redshift, while the second acts as a constraint on the speed of the spacecraft. It is important to stress that we 
have here a term associated with the energy flux, the last one of the left hand side in equation ( pi[ ) . 

Let us take as an example the case of an expanding wormhole, with Q = e 2uJt and domain for t in the interval 
(— 00, 00). Now, equations (p3|) and (pjh can be rewritten in the following compact form, 
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It is seen from these equations that in this expanding "wormhole universe" the tidal forces felt by an observer at 
fixed r diminish with time. The traversability criteria will be then satisfied for any value of oj if t is big enough. 
On the other hand, when the wormhole is contracting (i.e. fl(t) — e~ 2uJt ), the traversability criteria constrain the 
possible values of a, (3, and 4>tti more and more as the geometry evolves, eventually rendering the trip impossible. 
The behaviour of the tidal forces is then intimately related to the conformal factor f2(t), and one would generically 
expect that if there is an expansion (contraction), the tidal forces will diminish (grow) with t. 

V. FINAL COMMENTS 

We presented the complete analytical solution for a restricted class of evolving wormhole geometries by imposing 
a suitable form of the trace of the stress-energy tensor for matter. This, in turn, depends on the particular election 
of the conformal factor. With this solution we have analyzed different scenarios of WEC violation, extending in this 
way previous works by Kar and Kar and Sahdev 15 16]], in the sense that the redshift function in our case is not zero. 



It is worth noting that we have introduced a nonzero energy flux in the stress-energy tensor, which also renders our 
non-static solution more general than the ones presented before. 

We should also remark that one of the keys for the existence of these solutions is their non-flat asymptotic behaviour. 
In an asymptotically flat space-time, the existence of traversable wormhole solutions in GR is forbidden by the theorem 
of topological censorship j2l]] . We have also studied some traversability criteria concerning tidal forces that a traveler 
might feel. We found out that, since the geometry itself is evolving in time, there could be in some cases an expansion 
of the throat of the wormhole which may favor the diminishing of destructive forces upon the trip. 
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TABLE I. Functions involved in WEC inequalities 
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TABLE II. Possible conformal functions 



fi(t) T(t) /i f 2 f 3 Domain for t 
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a + l a A +t 2 (a 2 +t 2 ) 2 (a 2 +t 2 )' 2 (a A +t 2 ) 2 — CO < t < CO 
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FIG. 1. Behavior of the radial part of the first WEC inequality, gi(4>), near the wormhole throat. The curves are, from black 
to grey, the corresponding to a — 2, 5 and 10 and all them are presented for the case <f> th = 2. 
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Fig. 2 




FIG. 2. As in Fig. 1 but for the function gi (<j>), involved in the second WEC inequality. 



Fig. 3 




FIG. 3. As in Fig. 1 but for the function g3(<f>), involved in the third WEC inequality. 



8 



